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ABSTRACT 

We report a A > 1 cm search for rotational molecular absorption towards quasars, 
now possible with the upgraded Australia Telescope Compact Array (ATCA). The 
targets were PKS 0201+113, PKS 0336-017 and Q 0537-286, where known damped 
Lyman-alpha absorption systems (DLAs) could cause redshifted molecular absorption 
in the 12 mm band of the telescope. We place 3(t upper limits on any HCO’*' 0 —> 1 
absorption features of < 30 mJy per 3.4 km s~^ channel. The non-detections could 
be attributed to the inherent low metallicities in DLAs leading to generally low H2, 
and thus HCO'^, column densities. In general, the detection of molecular rotational 
transitions in DLAs could be further hindered by a lower than expected CO-to-H2 
conversion ratio, whether due either to photoionization of carbon or its relative under¬ 
abundance at high redshift. 

Key words: quasars: absorption lines-galaxies: ISM-radio continuum: galaxies- 
cosmology: early universe 


1 INTRODUCTION 

Molecular absorption lines at high redshift can provide an 
excellent probe of cosmological physics such as the cosmic 
microwave background, values of the fundamental constants 


Combes 

1996c 


1997 


2001 



arlv Universe (e.g 

. Wiklind & 

Drinkwater et al. 

1998 

). How- 


ever, such studies are limited to the 4 known high red- 
s hift molecular absorption s ystems, towards TX S 0218+357 
(Wiklind & Combes 1995), P KS 1413+135 ( Wiklind & 
Combes 1997), TXS 1504+377 ( Wiklind & Combes [1996c ) 
and PKS 1830-211 (Wiklind & Combes 1998). In the search 
for new systems, one systematic approach is to target high 
column density absorbers with known redshifts. A conve¬ 
nient sample is the damped Lyman-alpha absorbers, which 
have neutral hydrogen column densities Ahi > 10^° cm“^. 
In order to sel ect a sample, we pro duced a catalogue of all 
known DLAs (Curran et al 2002band shortlisted those 
which are illuminated by radio-loud quasars (i.e. those with 
a measured radio flux density > 0.1 Jy). From this sample 
of 57 we selected those which have 12 mm or 3 mm fluxes. 

Recently, we completed a search for molecular absorp¬ 
tion towards 11 DLAs at A ^ 3 mm with the SEST 15-m 
and Onsala 20-m telescopes which, apart from one tenta- 
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^ A version of t his r.ata.lop-ne is continua.llv undated on-li ne and 


f ive detection, only l ead to upper limits for 18 transitions 
( [Curran et al. 2002a). With the upgraded ATCA it is now 
possible to improve on these previous attempts. In this paper 
we present the results of our first search with this telescope 
- observations towards the known southern centimetre-loud 
quasars occulted by DLAs in which a commonly detected 
transition falls into the 12 mm band. 


2 OBSERVATIONS 

The observations were performed in June 2002 with the 
ATCA at Narrabri, Australia during excellent weather con¬ 
ditions which gave good phase stability. The telescope has 
recently been upgraded with the addition of 12 mm and 3 
mm receivers to antennae 2, 3 and 4 (see Wong & Melatos 
2002), thus permitting the search for redshifted molecular 


rotational lines with this instrument. As mentioned above, 
the 3 mm band has been used quite extensive l y in previous 


searches towards DLAs (Wiklind fc Combes 1994b, 


199E 


1996b; Curran et al. 2002aD, although to date no 12 mm 


is available from littp://www.phys.unsw.edu.au/~sjc/dla 


searches have been published. At this wavelength we are able 
to take advantage of the lower system temperature (~ 60 
K) and better atmospheric stability to observe towards the 
southern (<5 ^ 30°) quasars of sufficient flux^ occulted by 

^ Originally PKS 0336-017 was confused with Q 0336-019 (J 
0339—017), which was then used as a calibrator for the source 
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Table 1. The 22 GHz illuminated DLAs where the redshifted HCO”*" 0 —> 1 falls into either the 16.089-18.888 or 20.089- 
22.488 GHz ATGA bands, the DLA redshift and i/qUs is the expected frequency of the redshifted HCO't 0 —> 1 line, 

given to the number of significant figures available from the optical data. The red, visible and blue magnitudes of the quasar 
are given as well as the measured 22 GHz flux densities from both the literature (Lit.) and our observations (Obs.). Any 
difference between these two values may be due to quasar variability. 


Source 

Coordinates (J2000) 

^aJas 

^obs 

Magnitude 

Ri £ 

22 (Jy) 

Calibrator 

~ S 22 (Jy) 


h m s 

d ' " 



B 

V 

R 

Lit.Q 

1 Obs. 


Lit. 

Obs. 

0201+113 

02 03 46.7 

11 34 44 

3.38639 

20.331 

21.7 

19.5 

18.8 

0.55 

0.59 

0202+149 

2.00 

1.41 

0336-017 

03 39 00.9 

-01 33 18 

3.0619 

21.958 

19.5 

18.8 

19.1 

- 

0.15 

0336-019 

3.37 

2.2 

0537-286 

05 39 54.3 

-28 39 56 

2.976 

22.43 

19.8 

19.0 

18.9 

1.66 

0.58 

Self calibration 



The flux density for 0201+113 is from Terasranta et al 


f 2001) and for 0202+149, 0336-019 and 0537-286 from the ATGA 


calibrators site (http:/ /www.narrabri.atnf.csiro.au/calibrators/|). The uncertainty of ^ 15% in these calibrators is the major 


contributor to the errors in the measured flux densities. 


DLAs at redshifts of ~ 3 (Table In order to minimize 
the bias introduced by using such optically selected sources, 
ideally we would have selected the most visually faint (and 
hence dusty) objects. However, the ATGA still has a rela¬ 
tively narrow tuning range in the 12 mm band. We therefore 
selected all of the “high flux” DLAs which can be observed 
with this restriction (see Table |^. Fortunately, all three 
are relatively faint (by DLA standards) with 0201+113 and 
0537-286 also being red (B — R = 2.9 and 0.9, respectively). 

Because of the lack of precision in the optical redshifts 
(most significant for 0537-286, Table ^), we sacrificed one 
polarisation in order to overlap two of the widest available 
(64 MHz) bands. This enabled us to cover an uncertainty of 
« ±0.01 in redshift while retaining a relatively high spectral 
resolution of 3.4 km s”^. Finally, post observation, baseline 
2-4 was discarded because of phase referencing problems. 
No other flagging of bad data was required. For 0201 + 113 
and 0336-017 the bandpass of the calibrator was removed 
from the spectra and in the case of 0537-286, which was self 
calibrated, we removed a low order polynomial to “flatten” 
the bandpass. 


3 RESULTS AND DISCUSSION 

In Figs. 1^ to ^ we show the time averaged spectra over both 
good baselines and note that there are no HCO’*' 0^1 
absorption features of ^ 3a per 3.4 km s“^ channel in these 
DLAs. From the r.m.s. noise levels we derive optical depth 
limits for a resolution of 1 km s“^ and column density limits 
according to Nmm oc f rdv for an excitation temperature of 


10 K (see 


Curran et al 


2002a). These are listed in Table 


^ together with all previously published results. 

From the table we see that, after 3-5 hours per source, 
our limits are among the lowest of all the searches and 



Doppler—Corrected Frequency [GHz] 

Figure 1. HGO'*" 0 ^ 1 at z = 3.386 to.oig towards 
0201+113. The la r.m.s. noise is 10 mjy. The overlap be¬ 
tween the two spliced bands is 14 MHz. In this and Figs. ^ 
and ^ the arrow shows the expected absorption frequency 
according to the published DLA redshift (Table ^). 


we achieve the most sensitive search for HCO^ absorp¬ 
tion in a DLA published to date. However, there are un¬ 
certainties in the conversion ratio betweeen Njjqo+ 

Nco (and thus Ang which the CO traces^]). The ratio typ¬ 
ical of Galactic star form i ng cl ouds is Nco ~ 10"*Ajjqo+ 
(e.g. Wiklind & Combes 1995), whereas that for Galac¬ 
tic absorbers towar ds extragalacti c continuum source]^ is 
Nco ^ 10^Ajjco+ (Liszt & Lucat 1T98 ). For hig h redshift 
clouds, Nco > 10^ivHco+ dWiklind fc Combes 1995| ) or, 
more specifically, Nco = 7500Ajjco+ for one of the 4 known 
absorbers (Section 1 ), the gravitational lens PKS 1830-211 
(Menten et al. 1999). Therefore our results are of little use 
in assigning upper limits to the molecular hydrogen column 


(Table |]). The measured 22 GHz flux density of S 22 = 0.15 Jy 
now joins So. 4 = 1.31, Si. 4 = 0.60, S2.7 = 0.45 an d S.e; n = 0.30 


Jv for the measured radio flux densities of 0336—017 (Curran et al, 
2002b[). 

HCO+ is the strongest and most commonly detected molecule 
in the four known high redshift molecular absorption systems (see 
Wiklind &: Combes references). Note that no CO transitions fell 
into either of the two sub-bands at these redshifts. 


^ Although molecular hydrogen constitutes the bulk of the molec¬ 
ular gas in interstellar space, it cannot be observed directly due 
to its small dipole moment and moment of inertia. Since CO is 
the next most abundant molecule after H2, this “tracer” is ex¬ 
tremely useful in the study of the bulk ISM and is therefore the 
most studied molecule in external galaxies. 

^ This applies in the regime in which we are interested, i.e. CO 
column densities of Nqq ^ 10^^ cm“^ (Table ^). Above this the 
CO begins self-shielding and subsequently Nqq > 10^A’jjqq+. 
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Table 2. Summary of published searches for rotational molecular absorption in DLAs. Uohs is the approximate observed 
frequency (GHz), V is the visual magnitude of the background quasar, Ahi (cm~^) is the DLA column density from the 
Lyman-alpha line and r2i cm is the optical depth of the redshifted 21 cm Hi line (see Curran et al. 2002b). The optical 
depth of the relevant millimetre line is calculated from r = — ln(l — So’rms/'S'cont)? where Urms is the r.m.s. noise level at a 
given resolution and S'cont is the continuum flux density. This is done for a resolution of Au = 1 km s~^ (Tmm), where we 
quote only the best existing limit. Note that due to simultaneous flux measurements with the ATCA, unlike many of the 
other results given, we minimize errors due to variable fluxes (see Table For all optical depths, 3<j upper limits are quoted 
and ” designates where 3cr > S'cont, thus not giving a meaningful value for this limit. Blanks in the T2 i cm field signify 
that there are no published Hi absorption data for these DLAs. The penultimate column gives the best existing limit of the 
column density per unit line-width (not to be confused with Au) estimated for the transition [cm~^(km s“^)~^]. 


DLA 

■^abs 

Transition 

^obs 

V 

Nhi 

'^21 cm 

"^mm 

Nmm/dv 

Ref. 

0201+113 

3.38639 

HCO+ 0 

1 

20.3 

19.5 

2 X 10^1 

0.04 - 0.09 

< 0.1 

< 9 X 10“ 

7 

0235+1624 

0.52400 

CO 0 ^ 

■ 1 

75.6 

15.5 

4 X lO^i 

0.05 - 0.5 

< 0.06 

< 4 X IQi'^ 

2 


0.52398 

CO 1 ^ 

• 2 

151.3 




< 0.09 

< 6 X IQi® 

4 



HCO+ 3 

4 

234.1 




< 0.3 

< 4 X 10i2 

4 


0.523869 

CS 2 ^ 

3 

96.4 




< 0.9 

< 1 X IQi® 

6 

0248+430 

0.3939 

CS 2 ^ 

3 

105.4 

17.7 

4 X 10^1 

0.20 

- 

- 

6 

0336-017 

3.0619 

HCO+ 0 

1 

22.0 

18.8 

2 X 10^1 

< 0.005 

< 0.2 

< 2 X 10i2 

7 

0458-020 

2.0397 

HCO+ 2 

3 

88.0 

18.4 

5 X 10^1 

0.30 

< 2 

< 8 X 10i2 

5 


2.0399 



88.0 




< 0.4 

< 1 X 1012 

6 


2.0397 

CO 2 ^ 

■ 3 

113.8 




- 

- 

3 


2.0399 

CO 3 ^ 

■ 4 

151.7 




< 1 

< 2 X IQl® 

6 

0528-2505 

2.1408 

CO 2 ^ 

• 3 

110.1 

19.0 

4 X 10^° 

< 0.2 

No published 3 mm fluxes 

3 

0537-286 

2.974 

HCO+ 0 

1 

22.4 

19.0 

2 X lO^o 


< 0.08 

< 7 X IQii 

7 

0738+313 

0.2212 

CO 0 ^ 

■ 1 

94.4 

16.1 

2 X 10^1 

0.07 

- 

< 4 X lOi® 

6 

0827+243 

0.5247 

CS 2 ^ 

3 

96.4 

17.3 

2 X 10^0 

0.007 

< 0.4 

< 1 X 1013 

6 


0.52476 

CO 0 ^ 

• 1 

75.6 




< 0.2 

< 1 X lOi® 

8 

08279+5255 

2.97364 

HCO+ 0 

1 

44.9 

15.2 

1 X lO^o 


No published 1 mm fluxes 

6 



CO 2 ^ 

• 3 

87.0 




No published 7 or 3 mm fluxes 

6 

0834-201 

1.715 

HCO+ 2 

3 

98.6 

18.5 

3 X 10^0 


< 0.4 

< 2 X 1012 

5 



HCO+ 3 

4 

131.4 




< 0.6 

< 8 X 10i2 

5 



CO 3 ^ 

• 4 

169.8 




- 

- 

6 

1017+1055 

2.380 

CS 2 ^ 

3 

43.5 

17.2 

8 X lOi'^ 


No published 1 mm fluxes 

6 



CO 2 ^ 

• 3 

102.3 




No published 3 mm fluxes 

6 

1215+333 

1.9984 

CO 2 ^ 

• 3 

115.3 

18.1 

1 X lO^i 


- 

- 

3 

1229-0207 

0.3950 

CO 0 ^ 

■ 1 

82.6 

16.8 

1 X lO^i 


< 1 

< 6 X lOi® 

6 



CO 1 ^ 

■ 2 

165.3 




- 

- 

6 


0.39498 

CO 1 ^ 

■ 2 

165.3 




- 

- 

4 

1328+307 

0.69215 

HCO+ 1 

2 

105.4 

17.3 

2 X 10^1 

0.11 

< 0.7 

< 2 X 10i2 

4 



CS 2 ^ 

3 

86.9 




< 2 

< 4 X 1013 

6 



CO 1 ^ 

• 2 

136.2 




< 1 

< 3 X lOi® 

4 



CO 2 ^ 

• 3 

204.4 




- 

- 

4 

1331+170 

1.7764 

CO 0 ^ 

■ 1 

41.5 

16.7 

3 X 10^1 

0.02 

- 

- 

1 


1.7755 



41.5 




- 

- 

1 

1451-375 

0.2761 

HCO+ 1 

2 

139.8 

16.7 

1 X lO^o 

< 0.006 

< 0.6 

< 2 X 1012 

6 



CO 0 ^ 

• 1 

90.3 




< 0.3 

< 2 X 1013 

6 

2136+141 

2.1346 

HCO+ 2 

3 

85.4 

18.9 

6 X lOi'^ 


< 0.3 

< 1 X 10l2 

5 



CO 2 ^ 

• 3 

110.3 




< 0.6 

< 4 X 1013 

5 



CO 3 ^ 

• 4 

147.1 




< 0.8 

< 2 X loi® 

5 



CO 5 ^ 

• 6 

220.6 




- 

- 

4 


Refere nces: fll [Pakahara et al. j 1984 ): (2) Fakahara et al 


(1995); (5) Wiklind & Combes (1996b); (6) Curran et al 


( 1987 ); (3) Wiklind S.c Combes (1994a); (4) Wiklind & Combet 
^ 2002a| ); (7) This paper; (8) A. Bolatto (private communication) 


density, . In order to bypass this uncertainty, we use the 
best Nco limit (from Table |^, to estimate an upper limit, 
though we are still forced to assign an expected line-width 
to a non-detected feature. For three of the four known high 
redshift absorbers this is ~ 20 km s~^ (the FWHM being 
7 km towards PKS 1413+135Q, [Wiklind &; Comte 


1994a 


1995 


1996b, 1998| ). Also, as discussed in [Curran et al 


See http://www.phys.unsw.edu.au/~sjc/dla-figl.ps.gz 


( ^002a ), where non-detections are concerned the spectral 


resolution of the data have to be taken into account since 
these determine the r.m.s. noise and thus the values of r 
(Table 1). 

Nevertheless, using the best available optical depth 
limit (CO 0^1 towards 0235-1-1624, Table ^ gives Nco ~ 
7 X 10^® cm~^ for a 20 km s“^ line detected at a resolution 
of 1 km s“^. This becomes Nco 3 x 10*^® cm“^ for the 
same (barely resolve dl line at 10 km s~^ reso lution. So for 
Nh 2 ~ 10'*A^co (e.g. Wiklind fc Combes 1995|) , and assum- 
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Doppler—Corrected Frequency [GHz] 

Figure 2. HCO+ 0 ^ 1 at 2 : = 3.062 lo'ojy towards 0336- 
017. The \a r.m.s. noise is 6 mjy. The overlap between the 
two spliced bands is 12 MHz. 



Doppler-Corrected Frequency [GHz] 

Figure 3. HCO+ 0 ^ 1 at 2 : = 2.976 lo'oJ” towards 0537- 
286. The Ict r.m.s. noise is 8 mJy. The gap in the spectrum 
occnrs since there was no bandpass calibration and the fitted 
polynomial is unreliable at these frequencies. 


ing that molecular absorption lines have a FWHM of 20 
km s~ , we can assert that A^H 2 ~ 1% Nm for the DLAs 
deeply searched for rotational molecular absorption^. 

For high redshift (2 > 1.8) sources, the ultra-violet H 2 
lines are redshifted to optical wavelengths and indeed H 2 


of ^ 1% ( 

Ge & Bechtold 

1997: Brianand & Petitiean 

1998 

Levshakov et al. 20001: Petitjean et al. 

200 c; Levshakov et al. 

2002; Ledoux et al, |2002). In the low metallicity environ- 


ments typical of DLAs (fZn/H] ~ —2 t o — 1 , i.e. ^ 0.01 to 


0.1 solar, e.g. Prochaska et al. 2001), Liszt (2002) argues 
that even if the gas is cool, the ambient ionization is suf¬ 
ficient to suppress the formation of H 2 . This may offer a 
potential explanation for our non-detections. 

The under-abundance of heavier elements introduces 
another uncertaint y: the applic a tion o f Gal actic CO-to-H 2 


conversion ratios. Black et al. 


( 1987| ) and Chaffee et al 


(1988) suggest that the photoionization of carbon would re¬ 
duce the ratio of A’h, to Nco to a tenth of Galactic values 


No te that the Nh^/Nqq ratio may exceed 10® (Liszt & Lucas 
1998) which could increase this value to iS 10%. 


and, fr om an upper limit of Nco f rom the A-X molecular 
bands, Brianand & Petitjean (1998) find Nco/Nm < 10“® 
in DLAs, which is indeed only 10% of the local ratio. In addi¬ 
tion, at earlier epochs of chemical enrichment we would ex¬ 
pect that there is simply a lack of raw materials to form sig¬ 
nificant amounts of tracer molecules (e.g. CO, HCO’^, HCN 
and CS) even in cold dark clouds where H 2 readily forms. A 
reliable conversion ratio could be obtained by searching for 
CO absorption towards the DLAs in which H 2 has been de¬ 
tected. However, the detection of H 2 biases towards optically 
bright sources and hence an under-abundance of dust. Fur¬ 
thermore, none of the quasars illuminating the 7 known H 2 
absorbing DLAs have appreciable 12 or 3 mm fluxes (Curran 
et ah, in preparation). 

In summary, not only do the low metallicities reduce 
the molecular hydrogen content in DLAs, but they may also 
further lower the column densities of other molecular trac¬ 
ers, making the detection of these much more difficult than 
expected. It is clear that further statistics are required in or¬ 
der to identify the factors most crucial (e.g. column density, 
visual magnitude or metallicity) in determining the molec¬ 
ular content of DLAs. It would be of interest to compare 
the metallicities of the four known molecular absor bers, es¬ 
pecially in the highest redshift case (2abs = 0.886, Wiklind 


& Combes 199f). However, these lie along the lines-of-sight 


to quasars which, although bright in the millimetre regime, 
are too optically dim (V> 20 in all cases) to estimate the 
metal-to-hydrogen ratio. 
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